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ABSTRACT 

Kinematic precision is affected by errors which are the result of either 
intentional adjustments or accidental defects in manufacturing and assembly of 
gear trains. This paper explains a general method for the determination of 
kinematic precision of gear trains. The general method is based on the exact 
kinematic relations for the contact point motions of the gear tooth surfaces 
under the influence of errors. An approximate method is also explained. 

Example applications of the general and approximate methods are demon- 
strated for gear trains consisting of involute (spur and helical) gears, 
circular-arc (Wi Idhaber-Novikov) gears, and spiral-bevel gears. Gear noise 
measurements from a helicopter transmission are presented and discussed with 
relation to the kinematic precision theory. 
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Transmission error is a measure of the kinematic precision of gear trains. 
Transmission error is defined as the departure of the meshed gear pair (or en- 
tire gear train) from a constant ratio of angular motions as defined by the 
ratio of tooth numbers. It is true that in a gear pair each gear has a whole 
number of teeth and this defines the nominal ratio of angular positions between 
the two. But the instantaneous ratio during a meshing cycle can vary slightly 
from the nominal ratio. Transmission error is the measure of instantaneous 
variation from the ideal nominal value. The precision of gears was investi- 
gated by Litvin [1]*, Litvin and Gutman [2], anu Michalec [3]. 

When the mating teeth in a. gear train have profiles that transmit motion 
having no error, they are said to be conjugate pairs. In theory, it is pos- 
sible to select an arbitrary shape for a driving tooth and then to find a pro- 
file for the driven tooth which will give conjugate action, involute gear 
teeth happen to have the same form for driving and driven member teeth. A 
benefit of the involute form is that small errors in center distance between 
gears will not produce transmission errors. This is not true for other tooth 
profile forms [4, 5j. 

In general, there are many causes for transmission error, and they cannot 
be avoided in practice. Such things as shaft misalignment, profile error, 
tooth deflections under load, mounting location errors, and gear support 
deflections may combine to cause transmission error. The effects of trans- 
mission error are most often harmful. These are high vibration and noise, 
pitting and scoring of gear teeth, and reduced reliability of the gear train. 
Sometimes it is beneficial to introduce small intentional errors into the gear 


♦Numbers in square brackets denote references. 
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tooth profile to compensate for the probable occurrence of accidental and una- 
voidable errors in the assembled and operating gear train. Tip relief to re- 
duce dynamic loading and combined mismatch in spiral-bevel gears to reduce mis- 
alignment sensitivity are two examples of intentional errors which are bene- 
ficial. 

I 

Baxter has studied the effect of various types of misalignment on tooth j 

i 

contact in bevel and hypoid gears [6]. Townsend, Coy, and Hatvani have 

i 

examined gear train noise as a test of its precision during an intentional 

i 

loss of lubricant destruction test [7]. 

Of course, the effect of all errors (intentional or otherwise) on 
transmission error must be predicted analytically if the gear design process 
is to remain rational and not collapse into a confusing heap of empiricism. 

Errors of manufacturing and assemblage of gears induce kinematic errors 
in gear-drives which may be presented by the following function: 

( 1 ) 

Here is the angle of rotation of the driving gear, 1, 

aQ = (acij^, aci2 « • • •) (2) 

is the vector of errors, and 

a »2 = ^2 ■ ^2 

is the kinematic error of the gear drive, represented as the difference between the 
theoretical and actual angles of rotation of driven gear 2. 

In this paper two methods to determine Function (1) are presented: (a) a 

numerical method for computer solution and (b) an approximate geometric solution 
which leads to simple, accurate results in an analytical form. 
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NOMENCLATURE 

shortest distance between gear axes of rotation 
eccentricity vector of gear i 
axial settings of gears (see Fig. 2) 
angular velocity ratio, gear i, gear j 
number of teeth on gear i 

vector function representing unit normals of surface of gear i 
in fixed coordinate system (f) 

change in unit normal vector due to errors in gear i 

change in unit normal vector due to point motion relative to 
gear i 

change in unit normal vector due to transfer point motion with 
gear i 

absolute velocity of tip of unit normal vector of surface i 

similar to dni^^ and dni]} but velocities rather than 
^r _,tr 

displacements 

center of rotation of gear i 

geometrical center of gear i 

pitch point 

vector of errors 

components of vector of errors 

vector from origin to axis of gear rotation 

base circle radius of gear i 

vector function representing surface of gear i in fixed 
coordinate system (f) 
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( 1 ) 


X 


9i 

0 

»2 


velocity of contact point 
coordinate system i 

coordinate system rigidly connected to frame 
absolute displacement of contact point of gear i 

displacement of contact point due to errors in gear i 

displacement of contact point relative to gear i 

displacement of contact point in transfer motion with gear i 
tangent plane 

surface coordinate of gear i surface 
similar to ds|.^\ velocities 

rather than displacements 


angular position of eccentricity vector of gear i 

initial angular position of eccentricity vector of gear i 

rotation vector representing position change of gear axis of 
rotation 

surface coordinate of gear surface i 

change in kinematic error function as measured on shaft of 
gear i 

proportionality factor 
radius vector 
surface i 

angle of rotation of gear 1 
actual angle of rotation of gear 2 

theoretical angle of rotation of gear 2 
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kinematic error function 

pressure angle 

angular velocity of gear i 
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a»2 

T_ 


THEORY AND EXACT SOLUTION METHOD FOR KINEMATIC PRECISION 
In the process of motion the tooth surfaces of two gears, and I 2 
(Fig. 1), are in tangency if the following equations are satisfied; 







•»2^ 


(4) 


n^ ( Uj^ ,9 2^ ) a n^ (U2*^2*^2^ 

Here r|’^ is the position vector of the contact point on gear i; n|'^ is 
the surface unit normal vector at the contact point M; u., and «. are 
the surface coordinates of the gear surfaces; and is the angle of 
rotation of gear i. Subscript f denotes a coordinate system which is rigidly 
connected to the frame. 

For a gearset with kinematic errors, represented by aQ^ and AQ 2 , con- 
ditions for tangency may be expressed as 


r^' ^ ( UjjOj j jAQj ) * r^ ^ (u2»®2*^2*^^2^ 

n^ ^ ( Uj 2 ,aQj ) = n^ ^(u2«e2«92>^Q2^ 
Equations (6) and (7) yield the functions 

’2^'^r^^r^^2^ =f2(»i) Aip2(»2.AQ2,AQ2) 

u ^ 2 » ♦ ^Q 2 ^ ® ^ (v 2 * j ♦ ^Q 2 ) 


( 6 ) 

(7) 

( 8 ) 
(9) 
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The functions 

' (u^yO>) • u • (<f y t &Q 2 } • ® .j { 9 j , AQ 2 ) (i ■ 1*2) (10) 

represent the path of the contact point on gear surface corresponding to the 
meshing of gears with errors of manufacturing and assembly. Functions 

represent the path of the contact point on gear surface corresponding to 
meshing without errors. Comparison of functions (1C) and (11) yields the change 
of the contact point path induced by errors. 

Consider the solution of equations (4) and (5) and ( 6 ) and (7). Vector 
equations (4) and (5) yield only five independent scalar equations since 
n(l)| = = 1. These equations may be presented as 

^ 1 ^ 1 * ^2 2 ^ “0 (j*lf2>...»5) (12) 

It is assumed that 

The symbol C^ denotes that functions f . have continuous partial deriva- 

J 

tives of the first order (at least) by all arguments. 

It is assumed that equation system (12) is satisfied by a set of para- 
meters 


p(l) 


iU][ 




,9 


( 1 ) ) 
2 ^ 


(13) 


and that surfaces Ej anti E2 ’f’ tangency at point Surfaces 

and ^2 in point contact in the neighborhood of if by the set of 

parameters the following Jacobian differs from zero: 
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0(f jtt2.f3.f4.f5) 

0(Uj.«j,U2.«2»»2^ 


3fj afj 3fj 8fj afj 

Tu^ TuJ 15^ 

• • • • • 

• • • • • 

^^5 ^^5 ^^5 ^^5 ^^5 

3Uj 3«j 3 U 2 3«2 ^2 


If inequality (14) is satisfied, equation (12) may be solved in the 
neighborhood of with the functions 


( 14 ) 



j).e^(»j^).U2(»^).e2(»^).» 



Cl 


(15) 


The function ''“epresents the ideal law of motion. In most cases 

(for conjugate tooth action) function linear. 

Equations (6) and (7) also yield a system of five independent equations in 
six unknowns (uj,ej,<pj,U2.«2»*P2^ 

g j ( Ui j^i ,9 1 . U2.^2 2 * ^ (j “ 1.2, ....5) (15) 


It is assumed that this system is satisfied by a set of parameters 


p'2) 


(17) 


with the same value of as in the set If in the neighborhood of 


(21 

P' ' the Jacobian 


0 ( 91 . 92 . 93 . 94 . 95 ) 
0 ( ^1 .^1 . ^2 *®2*^ 2 ^ 


^ 0 


(18) 


then system (16) may be solved with the functions 


|u^(»^aQ),ej^((p^.aQ),U2(<y^,aQ),e2(?j^.aQ),(»2(»j^.aQ)| e qI ( 19) 
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Function represents the actual law of trotion transformation - 

the law of transformation of motion which corresponds to errors of manufac- 
turing and assembly. Kinematic errors of the gear drive are represented by 
the function 

a»2 " (20) 

This method of solution can provide, not only the kinematic errors of a gear- 
set, but also the new path of the contact point (see functions (10)). 

In general, the numerical solution of a system of five nonlinear equations 
is a difficult problem which requires many iterations. To save coin)uter time 
an effective method of solution was recently proposed by Litvin and Gutman 
[2]. The principle of this method follows: 

The system of equation (16) may be represented as 



» 0 

(21) 


» 0 

(22) 


= 0 

(23) 

^ 4 ! u 1 » ©i f U 2 *^2 2 ’ ^ ^ 


(24) 

^5^ *^1 * ^2 *^2 *^2 * ^ “ 0 


(25) 


Equations (21) to (23) are determined from vector equation (6), and 
equations (24) to (25) from vector equation (7). Here, A represents the 
shortest distance between the axes of rotation of the two gears and ana 
H 2 represent the axial settings of the gears (Fig. 2). 
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Systems Sj(Xj,Yj,Zj) and S2(X2,Y2tZ2) shown in figure 2 are rigidly 
connected to the driving and driven gears, respectively. Now suppose that some 
points and M2(u2,«2) surfaces and I2 chosen. 

With a set of known parameters (uj,©j^,U2.e2) * equations (24) and (25) become 
a system of two equations in two unknowns which may be expressed as 



(26) 

^2^^1*’2^ ■ 0 

(27) 


Upon solving for and ©2» checks that the following equations are 
satisfied; 


A — 1 2 * ^ ^ 

(28) 


(29) 

^2 ~ 2 ♦ ^2 ’^2 2 * ^ ^ 

(30) 


wijere A, Hj, and aQ are given values. 

In generil, the solution of the above two systems of equations ((26) and 
(27) and (28) and (30)) requires an iterative procedure. In practice, one of 
the four variable parameters (u^,©pU2,«2) fixed, and the other 
three are changed such that the two equation systems are satisfied. 

The advantage of the above method lies in the ability to divide the system 
of five equations ((25) to (30)) into two subsystems of two and three 
equations, and to solve them separately. 

This method was applied to investigate the sensitivity of Wildhaber- 
Novikov gears to errors of center distance mounting [1]. These gears are 
generated by two rack cutters which have normal sections as shown in Fig. 3. 
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Surfaces of the rack cutters are in tangency along a straight line N - M* 
which Is parallel to axis and passes through point M of the normal 
section. In the normal section, the shape of each rack cutter Is a circular 
arc of radius (1 - I, II). The location of point H Is defined by the 

parameter e • (Fig. 3). The line M - M' generates a helix on the gear 

tooth surface which Is the path of contact points. Although the procedure 
described above Is primarily for numerical solution. In this case analytical 
results were obtained. 

The Investigation showed that the error of center distance sA resulted 
In the change of location of the contact point path. The new location of the 
path Is represented by the equation 


Rj -y (pj sin e * r^)^ + (pj cos e sin (31) 

Where Rj Is the radius of gear 1 cylinder on which lies the helix of the 
new contact point path. Parameters e, pp Pp, and bp are shown In 
Fig. 3; r^ Is the pitch cylinder radius. 

With aA - 0, 


sin 



(32) 


where is the parameter corresponding to the desired location of the 
contact point path. 
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KINEMATIC RELATIONS BETWEEN PARAMETERS OF CONTACT POINT MOTIONS 
The following relations are the basis of the second method for the 
determination of kinematic errors In gear trains. As stated above, the 
tangency of gear tooth surfaces Is represented by equations (4) and (b). 
Because of the continuity of tangency of these surfaces, it Is required that 

^ *♦ j ^( u2»^2*^2^ (34) 

Here r^^^ (1-1,2) Is the velocity of the surface contact point In abso- 
lute motion (with respect to the frame); n^’^ Is the velocity of the tip of 
the surface unit normal In absolute motion (also with respect to the frame). 
Henceforth, r|’^ is designated as v^^^ and n|’^ as 

The velocity of absolute motion may be represented as the sum of two 
components: (a) the velocity of transfer motion (together with the surface) 

and (b) the velocity of relative motion (with respect to the surface). 
Consequently, 


,(1) . M) , ,(1). 

_abs ^tr ^r * 

v(2) . 

“abs 

v(2) 


(35) 

Habs 2tr 2r • 

Tabs “ 

^tr 


(36) 


For a surface represented by a vector-function 


OWQINAL PAGE I3 
OF POOR OUALITY 


and a surface unit normal 

,0) 




(i = 1.2) 

the following comes from definitions (a) and (b) above: 


(38) 


v(^) 

3r^^^ dv- 

9 

y(i) 

„(i) 

du^ 3r^^^ 

do- 

!tr 

“ Ht“ 


‘"i 

It 

30i 

It" 


an^^^ df. 



du^ 3n^’^ 

dOi 

» » j. 


3U. 

dr 

30i 

It“ 


(39) 


(40) 


Transfer velocity may also be determined in a kinematical way by supposing 

that a gear with surface rotates about an axis that does not pass 

through the origin 0. of coordinate system S^, which is rigidly connected 

to the frame. Vector is the vector of angular velocity of the gear's 

rotation. Then, 


^ X (41) 


where r^^^ is the position vector drawn from origin 0^ to the contact 
point on the tooth surface and is a vector drawn from 0^ to an ar- 
bitrary point on the gear's axis of rotation. 

The transfer velocity n^^ is represented by the equation 


n{’^ = X n^’^ (42) 


Equations (33) to (36) yield the following kinematic relations for two 
tooth surfaces which are in continuous tangency: 


(1) * (1) (2) * (2) 

V V = * V 


( 43 ) 
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= 


( 2 ) + ;( 2 ) 




Hr 


(44) 


Equations (43) and (44) were first proposed by Litvin [1], On the basis of 
these equations, important problems in the theory of gearing were solved, such 
as avoiding tooth undercutting, deriving the relations between curvatures of 
two gear tooth surfaces in mesh, and determining the kinematic errors of gear 
drives which are caused by errors of manufacturing and assembly. 

APPROXIMATE METHOD OF CALCULATION OF GEAR DRIVE KINEMATIC ERRORS 
As a general rule, kinematic errors of a gear drive determined by the 
exact method must be obtained in a numerical way using a computer. This is a 
disadvantage of the exact method. Therefore, an approximate method with the 
opportunity to obtain accurate results analytically is now presented. 

Figure 4 shows two gear surfaces and ^2 which are not in tangen- 
cy due to errors of manufacturing and assembly. Points M^^^ and M^^^ do not 
coincide, position vectors rj^^^ and rj.^^ are not equal, and surface unit 

normal vectors n^^^ and n^^^ do not coincide. To bring the two surfaces 
into contact it is sufficient to hold one gear fixed and rotate the other gear 
by an additional small angle. Since the gear with surface is the driv- 
ing gear, it is preferable to fix the position of surface Zj^ and rotate 
surface ^^2 to bring it back into contact with z^ The additional 
angle of rotation a» 2 represents the change of the theoretical angle of 

rotation 9 ^ which is exerted by errors of the manufacturing of manufacturing 
and assembly. The A 92 is as yet an unknown function of the vector of errors 
aQ and varies in the process of motion. Thus 

i^2 = f(i'i»AQ) (45) 
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The determination of the Function (45) is based on the following kinematic 
relations, which are analogous to (43) and (44): 


»4r’ ^ <'4” * “4^* 


, 42 ) , , 42 ), ds<2) 

_tr y y 


(46) 


dn{^^+ dn^^^ + dnj[^^ . dnJJ^ + dnj^^ + dn^^^ 

«tr ,r .q „tr 


(47) 


where ds^^^ is the displacement of contact point and dn^^^ (i « 1.2) is 
the change in direction of the surface unit normal due to errors of manufac- 
turing and assembly. To bring the surfaces into contact, it is sufficient 

to rotate only gear 2, holding gear 1 at rest. Therefore, and dn^j;^ 

are zero, and 




ds 


(2) . ds(2) 


tr 


ds'^> 

-.q 


(48) 


dn 


( 1 ) 




= dnf2^ + dn 
tr _ 


( 2 ) 




(49) 


To determine relations between ds|2)^ 
ing scalar products: 


and dSq2) 


take the follow- 


(ds 


( 1 ) 


* *^!q 


( 1 ) 


) = n . (ds 


( 2 ) 

tr 


+ ds 


( 2 ) 


ds(2)) 


(50) 


Since vectors ds^^^ and ds^^^ must lie the common tangent plane T, equation 
Jf 

(50) is reduced to the following: 


n . ds[J^ = n . (dsj|^^ - dsj^)) 


( 51 ) 


The vector ds|2^ may be represented by the following cross product: 
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(52) 


ds 


( 2 ) 

tr 




where is the incremental angle of rotation of gear 2 and r^^^ is the 
position vector drawn from an arbitrary point of the axis of rotation to the 
contact point M. 

Equations (51) and (52) yield 



Equation (53) is the basic equation for the determination of kinematic errors 
of gear drives. Its application will be demonstrated in the following sections. 

Analogous scalar products can be composed on the basis of equation (49). 

It can be proven that these scalar products are zero because the vectors in 
equation (49) all belong to the tangent plane. Hereinafter, the following 
notations will apply; 






(54) 


where zaq^-^^ and raq^^Vepresent the sum of linear-error vectors due to 
manufacture and assembly of gears 1 and 2, respectively. 

In many cases, however, errors in gear trains do not result from linear 
displacements, but rather from angular displacements. For instance, kine- 
matic errors may result from the misalignment of gear shafts. 

Figure 5 shows the axis of gear 2 rotation a-a in its ideal position. 
Suppose that, due to an error of assembly, axis a-a is rotated about a 
nonintersecting axis B-B. Such an error of assembly may be represented by 
the vector A6, which is directed along axis B-B, where the direction of 
A6 corresponds to the direction of rotation by the right-hand rule. 


of contact point M may 


With the given vector the displacement 
be determined as follows: 

(a) Vector aa, directed along the axis B-B, is replaced by an equal vector 
aa, which passes through the origin Op and the vector-moment R x aa. Here 

R is a position vector drawn from Op to an arbitrary point on the lire of 
action of vector aa (Fig. 5). 

(b) The displacement aq^^^ corresponding to aa may be represented by 

aq^^^ ■ aa X rp^^^ ♦ R x aa « aa x (rp*^^ - R) (55) 

A similar equation may be developed to determine the displacement of the contact 
point M exerted by an angular error corresponding to gear 1. 

With notations (54) the equation (53) for the determination of kinematic 
errors may be represented as follows: 

(av^^^ x r^”^ + raq) . n^”^ = 0 (56) 

where laq » “ £Aq|^^ and n^^^ is the surface unit normal at the contact 

point M. 

The location of the contact point M and the direction of the unit normal 
/Ml 

n' ' change in the process of motion. A further simplification of equation 
(56) results by assuming that in all positions the contacting tooth surfaces 
have a common normal which passes through the contact point M and the pitch 
point P. This is the fundamental law for uniform motion transmission. For 
planar gears the pitch point P coincides with the point of tangency of the 
pitch circles (gear centrodes). The pitch point for bevel gears is located on 
the line of tangency of pitch cones. In both cases the surface unit normal n 
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is col linear to PH and 

a + PH s + xn^^^ (57) 

Equations (56) and (57) yield 


X r^**^ + laq) . n^^^ . 0 (58) 


because 


(»<^l « r(">) . n<«) . 


4,(2) X (4'’) » xnW) 


= («‘2> X rj")) . n(") 


r.W) 




4«))] 


(.4^1 X rf> ) . n<")) 


(59) 


Application of equation (58) in place of equation (56) has the advantage 

(P ) 

that the location of the pitch point may be considered as a constant (r^ ' 

® const). However, the direction of the surface unit normal is a function of 
Three types of gears - involute (spur and helical) and Wildhaber- 
Novikov - are exceptions to this statement. For these gears the unit normal 
of the gear surfaces at their contact point does not change its direction. 

Because of kinematic errors, the angular velocity ratio fluctuates as the 
gear teeth pass through mesh. Figure 6 shows functions for two types of kine~ 
matic errors. The first is a piecewise, nonlinear, periodic function which 
has a period that depends on the ratio 


mi 2 


^1 


b 

a 


(60) 


where N- (i = 1,2) is the number of gear teeth and b and a are the 
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minimum integral numbers with which the ratio mj^2 can be expressed. The 

( 2 ) 

angle of rotation of gear 1 corresponding to the period of function Af' 
is equal to 2wa. 

Such functions of kinematic errors are caused by (a) the eccentricity of gears and 
(b) the crossing of the theoretical axis with the axis of rotation of gears. (The 
shortest distance between these axes rotates in the process of motion.) 

The second type of function ab^^^(bj^), shown in Fig. 6(b), has a period of 

of b, = This function is exerted by (a) errors in the generating process 
1 

of gear teeth and (b) errors of gear axis location which do not change in the 
process of meshing, etc. 


APPLICATION OF THEORY TO ECCENTRICITY OF INVOLUTE SPUR GEARS 


7 shows base circles of radii involute 

spur gears. The rotation centers of the gears are denoted 0^^^ and O 2 . 

If the centers of base circles 0^ coincide with centers of roation 0^ ^ 

(i = 1,2), and then vectors of gear eccentricity ae^ =0^0^ are zero. The 
involute curves are in tangency at a point M of the line of action KL. 

To model the meshing of gears with eccentricity, gears 1 and 2 are trans- 
lated from their theoretical positions by 0^ ^0^ = &e^ (i = 1,2). Now, the 
center 0^ will be offset from the center of rotation 0^^^. Because of 
this displacement of the gears, the tangency of their involute curves is brok- 
en: the curves will wither interfere with each other (intersect) or lose con- 

tact. To bring the involute curves into the contact once again, it is suffic- 
ient to rotate gear 2 by a small angle According to equation (58) 

the angle may be determined with the equation 


.( 2 ) 


(a 9 ^'> X r(P)) . n 


= (ae^ - ae2) . n^^^) 


The triple product results in (Fig. 7) 

(2) (P)^ JM) . (2)_(P) 

” l2 


(ap'"' X rV M . n' ' = A9 


2 cos *0 


(61) 


(62) 


ly 
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where is the pressure angle. 

Vectors of eccentricity ae^ and ae 2 form angles and »2 vector 

Q(i)Q(^). tj^ 0 se angles are measured in the direction of gear roatation (Fig. 7). 
The dot products yield 

aej^ . n^*^^ = ae^ sin(B^ + '{»q) 

ae^ . n^^^ = aeg sin(v^ - li^) (63) 

It results from equations (61) to (62) that 

,b(2) , “l * *0> * ' *o’ (64) 

- ;:p 


where 


r^^) = cos 


(65) 


is the radius of the base circle of gear 2. 

The center 0- (i = 1,2) of the base circle rotates in the process of meshing; 




(1) 


and O^^^are two instantaneous positions of this center (Fig. 8). Angles 


and B 2 can be represented as follows: 


®i “ ^10 ’r 


®2 " ®20 ’2 


(66) 


where Bj^q and B 2 Q correspond to the initial positions of centers 0^ and O 2 , 
with = (p^ = 0. 

Equations (64) and (66) yield 

^2j aej sin(9j + + ao2 sin(»2 ^ 2 ) 


av 


.( 2 ) 


(67) 


where 

n = (®10 ^ * 0 ^ ’ "^2 = (®20 - * 0 ^ 

For convenience, consider the kinematic error function to have zero magnitude at 
^(1) _ ,p(^) _ 0 . Then, the kinematic error becomes 
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Aej [sin(bj + tj) 

m— 


sin Yj] 


te^ [sinCb^ ♦ Y2) - sin Y2] 

mgi * jrj 


( 68 ) 


where 


m. 


'21 


“2 


“1 


r‘^> ' N 

*^b ”2 


Equation (68) represents the kinematic error of a gear train with two gears as 
the sum of two harmonics. The periods of these harmonics are equal to the 
periods of complete revolutions of the gears. 

Equation (68) may be made symmetric as follows: 


AO 


(2) ^ Aej[sin(9j + yj) - sin y^] Ae 2 [sin (»2 * ^ 2 ) ” sin Y 2 ] 


.( 1 ) 


m2i+ 


.( 2 ) 


"^2 


2 Ae^[sin(b. + y^) - sin y^3 


^ £ 
i=l 


.(i) 


m, 


'2i 


(69) 


Here, m 22 * 1 3nd *2 ~ *^21 

Equation (69) can be generalized for a train with n gears as follows: 


AO 


(n) 


J Ae.[sin(». + Y^) - sin y^] 

i=l Til 

^b 


m 


(70) 


where ao^^^ is the resulting kinematic error of the gear train represented 
as the angle of rotation of gear n (the output gear). 

A complicateu gear train is a combination of pairs of gears. The parameter 
Y^ may be represented as 
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for the driving gear of the pair, and as 
^i • ^0 - *0 

for the driven gear of the pair. For instance, for computational purposes, a 
train of three gears must be replaced by two pairs of gears. The idler (inter- 
mediate gear) is considered as the driven gear in the first pair, and as the 
driving gear in the second pair. 

Designate the kinematic error exerted by the eccentricity ae^ of gear 
number i as 


Ae.[sin(^. + Y.) - sin yJ 

; — '■ ( 71 ) 

where is the error of the rotation angle The maximum 

possible value of this error is 


AO- - AO. 

i,max i,min 


2 Ae. 


r 


(i) 

b 


The kinematic error of the train may be represented as 


(72) 


Ao^*^^ = r Ae. m . (73) 

i-1 ’ 

Usually gear trains are applied for the reduction of angular velocities 
and thus m^^ is less than 1. It results from equation (73) that the last 
gears of a train (numbers n, n - 1, n - 2) induce the largest part of the 
resulting kinematic error Therefore, the precision of these gears 

must be higher than the others. 
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Tlie largest value of the kinematic error function and its distri- 

bution above and below the abscissa depend on the combination of parameters 
T- (i - 1, 2, n). Figure 9 shows the distribution of a function ae^ 

(»^) exerted by eccentricity of gear i of the train. 

The resulting errors of a gear train may be compensated for in part, by 
definite rules of assembly of gears with eccentricity. For instance, for gears 
with tooth numbers N, » N 2 and equal eccentricities ae^ « ae 2 the 
resulting kinematic error will be approximately zero if eccentricity vectors 
ae^ and ae 2 (Fig. 7) are directed opposite each other. 

APPLICATION OF THEORY TO ECCENTRICITY OF SPIRAL BEVEL GEARS 
For spatial gears the word "eccentricity” is used to describe that the 
geometric axis of a gear is parallel to, but does not coincide with, its axis 
of rotation (Fig. 10). As the eccentric gear rotates its geometric axis gener- 
ates a cylindrical surface of radius ae. The eccentricity vector ae is a 
vector which rotates about the gear axis. The initial position of vector 
ae (its position at the beginning of motion) is given by angle a (Fig. 10). 

Figure 11 shows coordinate systems Sj(X^Yj,Zj) and S|r(X^,Y^,Z^), 
which are rigidly connected to gear 1 and the frame, respectively. System 
is an auxiliary coordinate system, which is also rigidly connected to 
the frame. Driving gear 1 rotates about axis Z^^. The position of ae^^ 
in coordinate system is given by the angle op whicli is made by 
ae^ aiid axis Xj^. The current position of ae^ in coordinate system 

(or $h) is defined by the angle (<p^ + o^) and the matrix equality 

. a,,] Ue‘>l] 

cos 0 sin Yj 

0 1 0 
-sin Yj 0 cos Yj 


ae^ cos(vj + o^) 
-ae^ sin(f|: ♦ a^) 
0 


(74) 
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Here Ue(»] is the matrix of vector ae^ in terms of its projections on 
on the axes of coordinate system The 3 by 3 matrix trans- 
forms elements of the column matrix Ae^^^ ^coordinate system 

from coordinate system S^. 

Matrix equality (74) yields 



] = 


Aej^COS(vj + aj)cOS 
-ae^ sin(e^ + o^) 
-aejCos(ej ♦ aj)sin 


(75) 


The vector of eccentricity of the driven gear can be defined in a similar 
way. Figure 12 shows coordinate systems and rigidly connected to 
gear 2 and the frame, respectively. The auxilliary coordinate system Sp is 
also rigidly connected to the frame. 


(21 

Vector ae^ ' is represented by matrix equalities 




( 2 )- 


cos Y- 


sin y. 


-sin Y' 


0 


0 cos Y- 


ae2 cos(»2 
ae2 sin(»2 ♦ 02) 
0 


(76) 


which after matrix multiplication gives 


[ae 



ae2 cos(f2 «2) cos Y 2 

A62 sin (»2 ♦ 02) 

aep cos(f2 * «2^ '*2 


(77) 


Kinematic errors induced by gear eccentricities may now be found by applying 
Equation (58) as follows: 
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‘^(2)p(2)n (2)1 . L\^) 



(78) 


»(1) *J2) 

(Fig. 12) is represented by the matrix 


Here aei*' and aei‘^ are given by matrices (75) and (77); vector afi 

-I 


( 2 ) 




] . [L,p] 


COS ^2 0 -sin ^2 

0 1 0 
sin Y 2 0 Y 2 

-a»2 sint2 
0 


0 

0 


(79) 


^2j 

Vector represents the position vector of the contact point which 

belongs to the line of action and n^ represents the common unit normal of the 
gear surfaces at their point of tangency. 

Equations (78) and (79) yield 


a»2 “ 




(80) 


-y cos ^2 * (x cos ^2 2 sin V 2 )'’y " ^ "^2 


where . ae^^^ - ae^^\ tae^ - ae^^^ - ae^^^ lae^ « ae^^^ - ae[^^ and 

the subscript f was dropped. 

Projections n^ and n^ of the surface unit normal n, and coordinates 
X, y, ano z of the contact point change in the process of motion. But since the 
changes in these variables are relatively small, they may be neglected (Fig. 13). 


n^ > sin c^. i^ ♦ cos c^ (cos B jf sin e k^) 
« 0, y^ « 0, z^ ■ L 


(81) 
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Equations (80) and (81) yield 

(«) 

^ L Sin >2 cos y cos a 

where 

A(»j) - a^ sin ♦ o^) + cos (?j ♦ oj) ♦ 

82 sin (f2 ♦ «2^ ^2 ^^2 

8j > -acj cos cos B 

*2 “ “^®2 *c ® 

bj • aej(cos yj sin i|»^ - sin yj cos sin b) 

b2 - -862(00$ Y2 sin ♦j. + sin Y2 cos sin b) 

»2 • jq 

It Is concluded from the form of equation (82) that kinematic errors In- 
duced by the eccentricity of spiral-bevel gears may be represented as the sum 
of four harmonics; the period of two harmonics coincides with the period of 
revolution of gear 1; the period of the other two coincides with the period of 
revolution of gear 2. 

The function 8*2(61) represented by equation (82) Is a smoothed, 
continuous function which serves as a first approximation. In reality the true 
function 62(61) breaks as different sets of teeth come Into mesh. This 
break can be discovered If 8*2(61) is determined by equation (80). 

GEAR TRAIN VIBRATION AND NOISE MEASUREMENT 

To illustrate the principles discussed on the subject of gear train preci- 
sion, Figs. 14 and lb are used. These figures show some frequency spectrum 
measurements made on a helicopter transmission running In a test stand [7]. 

The transmission had a spiral-bevel Input stage with IS teeth on the pinion 
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and 71 teeth on the gear. The pinion was turning at 6200 rpm and the output 
shaft at 3S5.5 rpm. The output stage was a spur planetary arrangement with a 
27-tooth sun, 3 planet gears, each with 35 teeth and a 99-tooth ring gear which 
was spllned to the transmission housing. An accelerometer was mounted on the 
case Immediately outside the spline. 

Figure 14 shows a broadband frequency spectrum measurement of the vibra- 
tion signal. The spur mesh frequency was 583 Hz and the spiral bevel mesh 
frequency was 1963 Hz. The spiral-bevel vibration signature was much stronger 
than the spur signature. This Indicates that the meshing accuracy was better 
for the spur mesh than for the spiral-bevel mesh. There are also other peaks 
in the spectrum at multiples of the fundamental frequencies of 1963 Hz ana 583 
Hz. These other peaks are the higher harmonics due to the noise and vibration 
pulsations as the teeth mesh being different from the pure sinusoidal shape as 
shown in Fig. 9. 

Figure 15, an expanded region of the autospectrum plot given in Fig. 14, 
shows many peaks that are symmetrically located about the spur gear mesh fun- 
damental frequency peak at 583 Hz. These peaks locate the sideband frequen- 
cies which are due to sources of modulation in the time dependent vibration 
waveform. Each source of modulation may produce one pair of sidebands If it 
is a harmonic modulator. If nonharmonic, the side bands will repeat many 
times, as in the case in Fig. 15. 

SUMMARY OF RESULTS 

Kinematic precision is affected by errors that are the result of either 
intentional adjustments or accidental defects in the manufacturing and assem- 
bly of gear trains. A general method for the determination of kinematic pre- 
cision of gear trains has been explained. The general method is based on the 
exact kinematic relations for the contact point motions of the gear tooth sur- 
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faces under the influence of errors. An approximate method was also explained. 

Example applications of the general and approximate methods were demonstra- 
ted for gear trains consisting of Involute (spur and helical) gears, circular 
arc (Wi Idhaber-Novikov) gears, and spiral-bevel gears. Gear noise measurements 
from a helicopter transmission were presented and discussed with relation to 
the kinematic precision theory. The following results were obtained: 

1. The exact numerical iterative procedure for finding kinematic errors, 

Ab 2 , is as follows: From equation (20) find 

where the angles *2 determined from an iterative 

solution of the nonlinear algebraic system of five equations, which are separ- 
able into two systems of two and three equations each as follows: 

Fi(»i, 92^ = 0, i = 1,2 
^ A — ( u , U2«02 2 ) ~ ® 

^ — K 2 ( u^ ,0^ ^ , U 2 , 02 ,^2 * ) ~ G 

^ H 2 - j ’ ^2 ’ ®2 ’^2 ’ ^ ~ ^ 

where A, H 2 , and aQ are given values and 92 ^5 determined by solving with 

aQ = 0, whereas is determined by the full solution of the five equations. 

2. The approximate equation for kinematic error is 


.(2) . .(P) 


(Ab' ' X ri ' X SAq) * n' = 0 


(M) 


where r^^^ is the radius vector to the pitch point, is the surface normal 

vector at the contact point M, and lAq is the sum of known error vectors. 
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3. Application of the formulas showed that W1 Idhaber-Novlkov gears are sen- 
sitive to any errors which cause changes In the center distance. A formula for 
the location of the gear- tooth contact point was given. 

4. It was found that for a pair of spur gears the kinematic error function 
due to eccentricities Is a sum of two single harmonics. For a multistage speed 
reducer. It was concluded that accuracy In the final stages has the most Impact 
on kinematic error. For gears with approximately equal known eccentricities* 
the kinematic error may be compensated for by directing the eccentricities 
opposite one another. 

5. For a pair of spiral-bevel gears, the kinematic error function due to 
eccentricities is a sum of four harmonics. 
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THE TOOTH DIRECTION 



ure 4 - Gear teeth not in contKt 






(a) Periodic with gear rotation, 
(b) Periodic with tooth mesh cycle. 
Figure 6. - Kinematic errors. 
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Figure 15. * Narrow band frequency spectrum showing sidebands 
around the spur mesh frequency. 


